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' Cosmic Microwave Background experiments must achieve very accurate calibration of their po- 

I larization reference frame to avoid biasing the cosmological parameters. In particular, a wrong or 

inaccurate calibration might mimic the presence of a gravitational wave background, or a signal 
from cosmological birefringence, a phenomenon characteristic of several non-standard, symmetry 
breaking theories of electrodynamics that allow for in vacuo rotation of the polarization direction 
of the photon. Noteworthly, several authors have claimed that the BOOMERanG 2003 (B2K) pub- 
lished polarized power spectra of the CMB may hint at cosmological birefringence. Such analyses, 
however, do not take into account the reported calibration uncertainties of the BOOMERanG focal 
plane. We develop a formalism to include this effect and apply it to the BOOMERanG dataset, 
finding a cosmological rotation angle a = —4.3° ± 4.1°. We also investigate the expected perfor- 
mances of future space borne experiment, finding that an overall miscalibration larger then 1° for 
Planck and 0.2o for EPIC, if not properly taken into account, will produce a bias on the constraints 
. on the cosmological parameters and could misleadingly suggest the presence of a GW background. 

PACS numbers: 98.80.Cq 
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^ ■ I. INTRODUCTION 

C/3 
. 

The latest measurements of the Cosmic Microwave Background (hereafter, CMB) anisotropies from satellite, ground 
■ based and balloon-borne experiments such as WMAP ([l|), B00MERanG-B2K (0]), ACBAR (0) and QUAD ([|) 
' have revolutionized the field of cosmology. The precise measurements of the acoustic oscillations in the angular power 
spectrum of the CMB have indeed not only provided a confirmation of the standard model of structure formation but 
I also constrained key parameters as the curvature, the energy density in baryons and cold dark matter, and the shape 
\^ . of the primordial inflationary perturbations. The new measurements expected by future experiments as Planck will 
T-H ' certainly improve our knowledge of the CMB anisotropies by providing maps at higher angular resolutions and in a 
(y-^ broader frequency spectrum (see e.g. Q). 

f*^ The highest impact for cosmology is however expected not from temperature but from polarization measurements. 

0^ , The current measurements of CMB polarization ([l|,0,[l|), while firmly establishing a detection broadly in agreement 
with the standard model predictions, lack the precision needed to improve substantially the cosmological constraints. 
^ Especially on large angular scales, an accurate future measurement of CMB polarization will certainly shed light 
on crucial aspects as the epoch and history of cosmic reionization ( 6]) and the amount of primordial isocurvature 
perturbations (7]). More ambitiously, the CMB polarization is perhaps the most promising tool for detecting the 
background of primordial gravitational waves, generated during inflation [^Q. 

The statistical properties of CMB linear polarization are indeed fully characterized by two sets of spin-2 multipole 
moments with opposite parities. In the standard inflationary paradigm the primordial magnetic-type modes {B 
or curl modes) are only produced by gravity waves (or "tensor" metric perturbations [lo|). while density (scalar) 
perturbations only excite electric-like parity modes {E or gradient modes) and they do not correlate with each other 
(see e.g. [ll[). A detection of polarization i3-modes at large angular scales will clearly hint towards the presence of a 
gravitational wave background and open a window on the physics of inflation, at extreme energy scales. 
However, since the seminal paper by Carroll et al. ((l^'l the possibility of a cosmological birefringence, or in vacuo 
rotation of the polarization direction of a photon, has been considered by several authors (see e.g. [Tsl . IT^ . [Tsl . [TgI . 

m, m, [13, HH and references therein). Even if exotic, cosmological birefringence is therefore predicted by a wide 
class of non-standard models and provides information on symmetry violations and on pseudoscalar fields beyond the 
standard model. The main effect of such rotation is to mix the E and B modes. Considering the usual C'i coefficients 
of the expansions in Legendre polynomials of the corresponding two-point angular correlation functions (see e.g. [llj]). 
the effect of cosmological birefringence is to correlate the spectra by an angle 9 such that: 
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Cf^ = CJ^cos(20)-CJ^sin(2^^) 

Cf^ = Cj^ sin(26i) + Cj^ cos(26') 

Cf ^ = Cf ^ cos2(20) + Cf ^ sin2(20) - Cf ^ sin(40) 

Cf ^ = Cf ^ cos2(20) + Cf ^ siii2(20) + sm{A0) 

cr ^ \ - en ^^nm + cr - .in'(2.)) (1) 



where the new quantities affected by birefringence are denoted with a prime. As we can see, even if one assumes 
a standard scenario where BB, TB and EB correlations are zero at the Last Scattering Surface (LSS), cosmic 
birefringence introduces a new BB component from EE. The above equations are exact for a constant rotation, and 
approximate very well the effects of a rotation angle varying with time if the time variation is sufficiently slow and 
the resulting total rotation after the propagation of photons (from LSS toward us) is sufficiently small, which is the 
case for most models of interest (see (l3l H. [isl [H, lH, [H, [13, IHI^)- The final result is that B- modes polarization 
could therefore not only be sourced by gravitational waves but also from a completely different physics. 

While both gravity waves and cosmic birefringence would represent sensational discoveries, one has to deal also 
with less exciting and definitely more common possibilities as experimental systematics. Cosmic birefringence could 
indeed be mostly exactly mimicked by a simple mismatch in the calibration of the principal axis orientation of the 
polarimeters. Let us remind that polarization can be fully described by the four Stokes parameters /, Q, U and V 
(see e.g. p3|). In particular, information about linear polarization^ are encoded in the Q and U parameters that 
depend on the reference frame. A wrong definition of the principal axis of the polarimeter which rotates the reference 
frame of an angle 9 will lead to un-propcrly defined and mixed Q and U, thus resulting in a mixing between E and 
B modes and producing a change in the polarization power spectra in the same way as described by Equations [TJ 

Cosmic birefringence and how to disentangle it from spurious systematics is therefore a very timely subject of 
investigation and has already motivated several recent papers. For example, the possibility of de-rotating CMB 
polarization and to reconstruct the rotation angle by measurements of higher-order TE, EE, EB, and TB correlations 
has been suggested in [2^1 and analyzed in detail in [25j, while in ^2^] the effects of beam systematics (which can 
mimic the rotation of polarization plane) have been considered. 

In this paper we focus on the effects that a rotation of the polarization plane, regardless of its origin, has on 
the determination of cosmological parameters from CMB data when the effect is not properly taken into account. 
In particular, a miscalibration can produce a bias on cosmological parameters as the baryon and cold dark matter 
densities tUb and Wc, the Hubble constant Hq, the optical depth r, the inflationary spectral index n^, and a false 
detection of gravitational waves other than obviously be confused with a genuine cosmological birefringence. We 
also discuss in detail the benefits of inserting a specific parameter describing a rotation of systematical origin in 
the analysis and we employ our formalism on recent data provided by the B00MERanG-B2K experiment as well 
on simulated datasets of the Planck and EPIC missions. In the case of B00MERanG-B2K we show that previous 
marginal detection of cosmological birefringence that have been pointed out by [27] , could be partially explained by 
a systematic effect. 

The paper is organized as follows: in section [IT] we discuss the calibration method adopted for the BOOMERanG- 
B2K experiment and few cases that can induce a systematic rotation of the plane of a polarimeter in future missions; in 
section lllll we set forth our methodology which we apply in section lTVl to revisit the B00MERanG-B2K constraints on 
cosmological birefringence; finally in section |V] we forecast the impact of a systematic rotation on future experiments. 



II. DETERMINING THE POLARIZATION ANGLE: SYSTEMATICS 



In general, polarimeters for the CMB must detect intensity and linear polarization. Both can be described through 
the Stokes I, Q and U parameters (see [1^). One horned radio-meters and Polarization Sensitive Bolometers (PSB) 
measure the quantity: 

Vt = y {/(7t) + eQm cos[2(<5 + (3)] + eU{%) sin[2(<5 + /?)]} + (2) 



^ As counterexample see llTl |22|| 

^ Circular polarization often neglected since not expected in the standard scenario 
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where t labels time, K is the detector gain, 7 the observed sky direction'^, n is the detector noise, 6 is the angle between 
the telescope reference frame and a meridian in the sky, and e and f3 are the polarimeter properties, so defined: 

• e is the cross polarization, the ratio between the minimum and the maximum signal from the polarimeter, when 
a linearly polarized source is rotated in front of it; 

• /3 is the angle at which the polarimeter is maximally sensitive to a vertical linear polarization, in the telescope 
reference frame. 

These parameters must be accurately calibrated. An incorrect measure of e will result in a fake attenuation or 
amplification of the polarization and therefore will induce a bias in the Q and U maps. An incorrect measure of /3 
will result in leakage from Q to U and from U to Q, thus to a further "mixing" of the CMB polarization angular 
power spectra in equations [1] above: any anaccounted rotation of the polarimeter will mimic a genuine rotation of 
the photon polarization plane. There are several techniques to measure the parameters f3 and e, depending on the 
instrument characteristics. 

In the case of B00MERanG-B2K, the polarimeters have been calibrated in the lab prior the flight. The calibration 
is obtained using a dedicated apparatus, able to illuminate the full telescope with a linearly polarized source, 
modulated by a chopper wheel. A description of this system is reported in [2^ (see Figure 12 in that paper). 
The polarization angle of the source can be rotated as needed to cover the full range (0 — 180)°, and the angle (3 
for each of the detectors has been measured by means of a best fit procedure. Even if the associated statistical 
error was small, there are systematic errors in the alignment of the calibration source and of the telescope during 
the calibration. These, however, affect in the same way all the 16 detectors of the instrument. So if we consider 
the dataset of the 16 measurements of /3 [1^, the average value of the difference between nominal and measured 
orientations is (/3) = —0.84° ± 0.69°. We assume that this is due to the above mentioned attitude misalignment. 
This value will be used as a prior for the /3 angle in the following analysis. Instead, the standard deviation of the 
differences between nominal and measured orientations is cr(/3) — 2.75° if we include all the 16 detectors, and is 
1.56° if we just include only the eight 145 GHz detectors. These values indicate the effect of rotation errors in the 
assembly of the polarimeters, and can be used as an estimate of the precision of the direction of the main axis for the 
B00MERanG-B2K polarimeters. The axis direction error propagates on the measured power spectra following eq.l. 
So in the following we will mainly consider the effect of a general rotation of the order of {(3) . 

In the case of satellite mission such as Planck calibration procedure are even more delicate. The optical characteristics 
can change in space due to thermo-elastic effects, and ground based measurements must be extrapolated to flight 
conditions. Alternatively these parameters can be measured in flight by observation of a strong linearly polarized 
source. The most interesting source from this point of view is the Crab nebula, which must be observed in advance 
by some ground based telescope in order to provide a reliable calibrator. In the case of Planck HFI, the issue of the 
polarimeter angle is complicated even more by the fact that, in the PSB employed, the polarimeter angle depends on 
the details of the lithographic properties of the detector and its coupling to the feeds. For the radio-meters on board 
Planck LFI e is expected to be very close to unity, given to the excellent isolation provided by the OMT. On the 
other hand, /3 is subject to similar potential uncertainties as bolometers and must be calibrated in flight. For Planck, 
the uncertainty on the polarimeter calibration properties are not been disclosed yet. In the following, we assess what 
calibration accuracy should be achieved with Planck to yield unbiased constraints on cosmological parameters. 



III. METHOD 

To understand the impact of a systematic rotation of the CMB linear polarization plane on the determination of 
cosmological parameters, we analyze a set of CMB data with a modified version of the publicly available Markov 
Chain Monte Carlo package cosmomc [sO*]. We consider 5 chains with a convergence diagnostics done through the 
Gelman and Rubin statistic, as now common in the literature. For the CMB data-sets we analyze either the already 
published B00MERanG-B2K data, including the TB and EB spectra from the NA data analysis pipeline (see 0]), 
either simulated mock data for the ongoing /future Planck 3 and EPIC[| satellite missions. For B00MERanG-B2K 
we consider two rotation angles (/3 and a) in order to take into account both a systematic and a possible additional 
angle of true, cosmic birefringence, origin. 



More accurately, the signal measured in a given direction 7 is given by a convolution between the true sky and an optical response 
function (or "beam"), which is a tensor in the case of polarization sensitive detectors. Such detail is irrelevant here, so we simplify the 
description assuming an infinite resolution scalar beam 
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For Planck and EPIC we focus on a single rotation angle and we study the impact on the cosmological inferences 
drawn by those missions and, in particular, on the detection of a i?-modes signal from primordial gravity waves. 

When simulating the Planck and EPIC future data we assume as fiducial model the maximum likelihood WMAP5 
parameters but with different values of miscalibration angles. We therefore sample the following eight-dimensional 
set of cosmological parameters, adopting flat priors on them: the physical baryon and Cold Dark Matter densities, 
LUb — flbh^ and luc — i^ch^, the ratio of the sound horizon to the angular diameter distance at decoupling, 9s, the scalar 
spectral index n^, the overall normalization of the spectrum As at k = 0.05 Mpc~^, the optical depth to reionization 
r, the tensor to scalar ratio of the primordial spectrum r and the systematic rotation angle /3. We also use a cosmic 
age tophat prior as lOGyr > to > 20Gyr. Furthermore, we consider purely adiabatic initial conditions, we impose 
flatness and we treat the dark energy component as a cosmological constant. 

In the case of the B00MERanG-B2K dataset we include an additional rotation angle but we fix r = 0.09 consistent 
with the WMAP5 measurement ([1]). 

IV. ANALYSIS OF THE BOOMERANG-B2K POWER SPECTRA 

Here we consider the TT, TE, EE, BB, TB, EB power spectra of the Cosmic Microwave Background obtained 
by the B00MERanG-B2K experiment flown in Jan. 2003 ([2])- This dataset has already being analyzed in previous 
papers. In particular, in [2^ . |16l | and [l9j a claim for a ~ 1.2cr indication for cosmic birefringence from this data was 
presented. Here we re-analyze the issue by including an additional systematic uncertainty f3 on the rotation angle 
together with the rotation angle a assumed as sourced by a true signal of cosmic birefringence. 

While both rotation angles affect the CMB spectra as in Eq. 1, there are few differences. First of all, only photons 
coming from the last scattering surface at redshift z ~ 1100 are affected by the birefringence angle a. We therefore 
do not rotate with a the CMB spectra for £ < 23 since the CMB signal from those angular scales is assumed to be 
sourced from smaller red-shift z ~ 10 when the reionization process starts (see discussion in jH). Second, we do not 
let the systematic angle (3 to vary as a free parameter but we impose a gaussian prior on it as = —0.9° ± 0.7° at la 
motivated by the discussion presented in sec. HIl 

Without systematic effects (/3 = 0) we found that the B00MERanG-B2K data, with a prior on the optical depth 
T = 0.09, yields a ~ Icr evidence for cosmological birefringence with a = —5.2° ± 4.0° at 68% c.l., in agreement with 
the previous results of [H, [l^ [13] ■ 

We have then included the possibility of a systematic angle (3. Taking into account /3 and the prior on it we instead 
find a = —4.3° ±4.1° at 68% c.l., marginalizing over /3. As we can see, even if the the weak indication for cosmological 
birefringence is still present, the inclusion of /3 reduces the statistical significance at about 1 — a level. 
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FIG. 1: Joint 2-dimensional posterior probability contour plot in the a-/3 plane, showing the 68% and 95% confidence level 
contours from the BOOMERanG-B2K data. 

It is also interesting to note that from Fig[Tl where we plot the 2-D likelihood contour plot from the BOOMERanG- 
B2K data on the a-/3 plane, the probability that the detected rotation is due only to a miscalibration (a = 0°,/5 = 
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—5.2°) is suppressed. As we can see in Table |T] the constraints on the cosmological parameters from BOOMERanG- 
B2K data are not much affected by the cahbration error, only the upper limit on r is slightly relaxed. 

Parameter Not including /3 Including j3 

7W 0.0228lHg^^ 0.02271^^8 

n.h^ 0.1231°°^^ 0.122iHf2 

Ho 74li lAtlt 

ns O.SSl"-^^! 0.88i»-^2? 

log[10io^,] 3.08l°:l« 3.08t°:i« 

r < 2.6 < 3.4 

-5.21^:^ -4.3^?;^ 

TABLE I: 95% confidence level on cosmological parameters obtained by analyzing BOOMERanG-B2K data with and without 
including the systematic rotation angle (3 in the analysis. 



V. FORECAST FOR PLANCK AND EPIC 



To forecast the impact of a systematic rotation of the CMB linear polarization plane on future CMB experiments 
we create mock datasets with noise properties consistent respectively with the Planck mission (see and the EPIC 
satellite proposal (see Q). For the simulated Planck datasets we consider the 100, 143, and 217GHz HFI detectors 
while for EPIC we use the single 150GHz channel. The experimental specifications we use are reported in Table HIl 
We consider for each channel a power noise = (Oa)^ where 9 is the FWHM of the beam assuming Gaussian profile 
and a is the sensitivity AT/T both from Table |TT1 We therefore add to each Ci spectra a noise spectrum given by: 



Ne^w-^expil{l + l)/ll), (3) 

where lb is given by h = v81n2/0. The assumed fiducial models are based on the WMAP-5 maximum likelihood 
parameters with no primordial gravitational waves background. However, wc also assume the possibility of 
a systematic rotation angle /3. For Planck we consider four fiducial models with f3 = —1°, —2°, —3°, —3.8°, while for 
EPIC, that will measure polarization with greater accuracy, we consider P = —0.1°, —0.3°, —0.5°, —0.7°. We choose 
those range of angles since, as we will see below, angles below ~ 1° will not affect the cosmic parameters in the case 
of Planck, while ~ 4° is the current constraint on cosmic birefringence (see [l[). In the same way, angles below ^ 0.1° 
will not affect the cosmic parameters in the case of EPIC, while ~ 1° will be presumably the conservative constraint 
on cosmic birefringence at the launch time of this satellite. 

For each of these datasets we perform two analysis: first we include the systematic rotation in the montecarlo 
allowing f3 to vary, then we redo the same analysis without including (3 among the parameters (i.e. fixing /3 to 
zero). In the first case we expect to recover the input fiducial model, in the second we study the possible bias on 
the cosmological parameters when the systematic miscalibration is not considered. We place particular attention on 
r, since an indication for r > would hint for a detection of gravitational waves, while no gravitational waves are 
assumed in the input fiducial model. 

Our results are reported in Table IIIII where we show the constraints on the standard cosmological parameters that 
could be derived from Planck HFI and EPIC respectively with and without considering in data analysis a possible 
systematic rotation of the polarization plane. 

As we can see in Table IIIH the constraints on cosmological parameters based on Planck mission data, when all 
the 3 HFI channels are considered , will be affected by mismatches when /? is greater than 1° ii (3 is not included 
among parameters in the analysis. Instead, as it is shown in Fig. [21 rotation angles of about or lower than 1° will 
have a minimal impact on the recovered parameters. This result on one hand is reassuring since this is the level of 
systematic error expected for this mission but on the other hand it is important to stress that, as will be explained 
below, there could be an intrinsic rotation of the power spectra due to cosmological birefringence, which could 
increase (decrease) the total rotation. However, larger values may bias the result towards lower values for the Hubble 
constant and baryon density and larger values for the cold dark matter component. Moreover, a miscalibration angle 
of about ~ 2.5° or larger could mimic a detection of B modes at more than 95% c.L. When the rotation angle is 
considered in the analysis the input cosmological parameters are correctly recovered. Finally, let us note that the HFI 
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Experiment Channel FWHM AT/T 





Planck 


70 


14' 


4.7 


fsky 


= 0.85 


100 


10' 


2.5 






143 


7.1' 


2.2 






217 


5' 


4.8 




EPIC 


150 


5' 


0.81 


fsky 


= 0.85 









TABLE II: Planck d and EPICjll experimental specifications. Channel frequency is given in GHz, FWHM in arcminutes and 
noise per pixel in 10~^ for the Stokes I parameter; the corresponding sensitivities for the Stokes Q and U parameters are related 
to this by a factor of \/2. 



Parameter 




Us 


r 


Ho 




r 


P 


Fiducial Model 


0.02341 


0.986 


0.090 


75.3 


0.1046 





See below 


Planck HFI 


Including (3 
















P = -1° 
P = -2° 
13 = -3° 
13 = -3.8° 


U.UZOOO_Q QQQ3Q 

n noqq7+o ooo33 

U.UZOO(_Q 00032 

n n9'<ss+o-ooo33 

U.UZOOO_0 00033 

n n9'?'?s+o ooo32 

U.UZOOO_g 00032 


n QCK9+0.0085 
u.yooz_0 g0g0 

0.9853lg:r8? 
0.985llg:»fo 
0.9852l»:r8^ 


0900+'^ °^^ 

0.090tHiS 
0.0895l°:°J}o 
0.0897l°:°099 


75.011:^2 
75.0tl:^ 
75.0tl:^ 
75.0^1:^ 


l049+° °°2* 

U.iUy:»_0 0025 

1049+°°°2'* 

U.lU'±y_0 0024 

1049+°°°2'^ 
u.iuy:y_0 0025 

1049+°°°^'^ 
u.iuy:y_0 0025 


< 0.022 

< 0.021 

< 0.022 

< 0.021 


-0 998+" °®^ 

U.»»O_0 og3 
Z.UU1_0 og0 

-3 002+° °*'' 

O.UUZ_0 ogy 

o sni +0.084 
O.OU1_0 032 



Not including 

f3 = -1° 
f3 = -2° 
(3 = -3° 
(3 = -3.8° 



+0.00032 
0.00032 
+0.00032 
0.00032 

02346+° °°°32 

U.UZOW_0 00032 
+0.00032 
0.00032 



0.02339 
0.02341 



0.02351 



0.9852l°:r8l 
0.985llg:rsl 
0.98561^^3 
0.9858to;oo83 



0.0898to:oo98 
0.0898to;oJg7 



75.0 
74.9 
74.7: 
74.6; 



+ 1.3 
1.3 
+ 1.3 
1.3 
1.3 
1.3 
1.3 
1.3 



1050+°°°25 

U.iUOU_0 0025 
in^'H+0 0025 
U.iUOO_0 0024 

in'S 7+0 0025 

U.iUO ( _n 0094 

0.1061 



0.0024 
0.0025 
0.0025 



< 0.026 

< 0.039 
0.034t°°l? 



EPIC 



Including P 
/? = -0.1° 
P = -0.3° 
P = -0.5° 
P = -0.7° 



0.02341 



+0.00015 
0.00015 

U.UZOW.o 00015 



0.02341: 



I +0.00015 
-0.00016 

02340+° °°°i'5 

U.UZOW_o 00015 



0.0045 
0.0048 

0.98571^^9 

+0.0047 
0.0046 
0.0048 
-0.0047 



0.9857: 



0.9859: 



-0.0067 
-0.0061 

n nsss+0-0067 
U.U888_0 ooeg 

-0.0068 
-0.0065 
-0.0068 
-0.0067 



0.0890; 



0.0891; 
0.0890; 



+0.51 
0.50 



75.26 
75.25tO;f3 
75.28«:f2 
75.23±°:a 



0.10444 
0.1044 
0.1044 
0.1044 



0.0010 
0.00098 
+0.0010 
0.0011 
+0.0011 
0.0011 
0.0011 
0.0010 



< 0.00072 

< 0.00069 

< 0.00071 

< 0.00070 



-0 0999+° °°**'' 

u.ua»a_o 00g7 

_n 9999+0 0092 
u.zyyy_o 0094 

-0 5000+° °°°'' 

U.OUUU_o 0091 

-0 6999+0 0095 
u.oyyy_o 0093 



Not including 
P = -0.1° 
P = -0.3° 
P = -0.5° 
P = -0.7° 



0.02340 



+0.00015 
0.00015 

n9S34+° °°°''5 
u.uzooy:_0 00015 



0.02324; 



+0.00015 
0.00015 

0.02308lHggl^ 



+0.0046 
0.0048 
n QS':!l+0 0049 
U.y»rfi_0 0049 

+0.0044 
0.0046 
n 0717+0.0048 
U.y I 1 I _0,oo48 



0.9855 



0.9784 



+0.0065 
-0.0063 



0. 

0.0895l°°°f6 
0.0908 
0.0925 



0.0066 
0.0067 
0.0063 
0.0075 
0.0067 



75 9ts+°-5i 
ro.zo_o 53 

74.76t°:g 

73.82«:1° 

79 /19+0.53 
IZ.^Z_o 51 



-0.0010 
-0.0010 

n 1054+°°°" 

U.iUO4_0 0011 
-0.0011 
-0.0010 
-0.0011 
-0.0011 



0.1045 



0.1074 
0.1104 



< 0.00074 

< 0.0011 
0.00086 
0.00129 



+0.00078 
0.00056 
+0.00091 
0.00068 



TABLE HI; 95% confidence level on cosmological parameters obtained by analyzing Planck HFI (100-143-217 GHz) and EPIC 
mock data with different systematic rotation angles. As we can see the input parameters are always recovered when P is 
included in the analysis. 



experiment will be able to measure a systematic angle when considered as an extra parameter, with a precision 
of cr(/3) ~ 0.1°. Since the expected calibration error will be of the order of ~ 1° it is clear that the use of a markov 
chain methods such those presented here will definitely be useful in understanding the experimental calibration. 
Let us now consider the case for the EPIC satellite proposal. As we can see in Table Iml and Fig. [3l the sensitivity 
of this experiment to P is higher respect to Planck, however, rotation angles of about or lower than 0.1° will have 
no effect in recovering the cosmological parameters. If we do not include /3 in the data analysis, larger values of the 
systematic angle will have a dramatic effect on all the cosmological parameters, moreover a miscalibration angle of 
about ~ 0.4° or larger could mimic a detection of B modes at more than 95% c.L. Finally the EPIC experiment will 
be able to measure a systematic angle /3, when considered as an extra parameter, with a precision of cr(/3) ~ 0.01°, 
an order of magnitude better than Planck HFI. 
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Channel Calibration angle /3 r (including P) r (without P) 



Planck 70GHz 


-1.0 


< 0.27 


< 0.27 


Planck 70GHz 


-3.8 


< 0.27 


< 0.28 


Planck lOOGHz 


-1.0 


< 0.082 


< 0.087 


Planck lOOGHz 


-3.8 


< 0.082 


< 0.13 


Planck 143GHz 


-1.0 


< 0.035 


< 0.039 


Planck 143GHz 


-3.8 


< 0.036 




Planck 217GHz 


-1.0 


< 0.077 


< 0.079 


Planck 217GHz 


-3.8 


< 0.078 


< 0.12 



TABLE IV: 95% confidence level Umits on r. 



We analyze in more detail the case of the Planck mission in Table ITVl where we show the effects on r separately for 
different channels. For the LFI 70 GHz channel, the results on r are not dramatically affected even for an unaccounted 
miscalibration oi f3 = —3.8°. The HFl 143 GHz channel has higher sensitivity and is expected to drive accurate results 
on cosmological parameters, even in the case where a set of frequency maps is used to provide a reference foreground 
cleaned CMB map. The 143 GHz channel can, however, be heavily affected by a miscalibration in the polarization 
angle. For instance, a value of /? = —1° will affect the expected bound on r by relaxing it at about ^ 10%. A larger 
miscalibration by f3 — —3.8°, if not properly accounted for, will have dramatic implications for this channel producing 
a fake detection of non-zero B modes at more than 2a confidence level. For the other two HFI channels considered 
here (100 and 217 GHz) there is no spurious detection of gravitational waves even for a systematic angle as big as 
—3.8°, even if the bound on r gets relaxed. Even a miscalibration f3 < 1° could, in the case of Planck, yield biased 
results for r if a genuine cosmological rotation (a) is present but not accounted for in the analysis. This because an 
unaccounted rotation can in principle bias r regardless of its origin (i.e. genuine or systematic). As a consequence, 
it is always necessary to include (3 among the derived parameters even in the case where the systematic rotation is 
kept under control. If a rotation is detected, disentangling its origin (cosmological or not) only based on the data is 
more troublesome. If cosmological birefringence is expected to scale with energy, as predicted by some models (see 
,[3 [3) the analysis of multifrequency data can shed light on the nature of the effect observed, since calibration 
errors are not normally expected to mimic this behavior. Models that do not display frequency scaling, on the other 
hand, cannot be constrained in this way. Another option is to rely on the i scale dependency of the rotation. Since 
the universe has undergone reionization, a fraction of the CMB photons have suffered an effective last scattering 
at z < 20, thus the magnitude of the accumulated rotations is expected to be severely suppressed at the lowest 
multipoles. A systematic rotation, instead, does not display the same behavior and hence this fact could potentially 
be used to disentangle the two effects. However, when we tried to verify this method on our Planck simulations we 
discovered that the number of multipoles that are affected by reionazation is too low for the method to be useful at 
disentangling the two contributions. Specifically, we chose not to rotate the multipoles ^ < 23, to find that the a — j3 
degeneracy can only be broken this way for amplitudes of order of several degrees, which is totally useless in a high 
precision context. Thus, with the possible exception of exotic models with small scale angular variations, models that 
do not exhibit frequency scaling can only be constrained assuming an accurate enough calibration of the telescope's 
focal plane. 



VI. CONCLUSIONS 



CMB polarization will be measured with exquisite precision by future experiments. A possible and worrying 
systematic that may affect this measurement is a miscalibration of the polarization reference frame. In this paper 
we have first investigated the impact of this systematic on the recent claim for a cosmic birefringence signal in the 
B00MERanG-B2K data. We have found that, when the calibration error in the polarization signal is properly taken 
into account, the constraint on a cosmic birefringence angle is a = —4.3° ± 4.1° at 68% c.l., i.e. almost consistent 
with a statistical fluctuation at one standard deviation. 

We have then analyzed the effect of a systematic rotation angle (3 on the determination of cosmological parameters 
from future CMB surveys, with a particular attention on the implications for a detection of a GW background. We 
have found that for the Planck mission, a miscalibration smaller than ~ 1° will have a negligible impact on the 
constraints for all parameters, while it could introduce a sizable bias for larger angles. For example, systematic error 
in the calibration of about ~ 3° (i.e. at about 3cr from the expected calibration of ^ 1°) would result in a ^ 3(7 GW's 
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FIG. 2: Likelihood distribution functions for Q.bh? , 0,ch^, Ho and r with (black solid line) and without (other solid lines) 
including the rotation angle (3 in the analysis. The noise considered for the mock data is equivalent to the overall expected 
sensitivity (3 channels) of the Planck mission. The rotation angles considered in the fiducial model are: /3 = —1 (red), [3 = —2 
(blue), [3 — ~3 (green) and /3 = —3.8 (cyan) 



detection. We have shown that when including in the usual MCMC analysis the additional parameter /3 to consider 
a possible rotation of the spectra (of systematic or cosmological origin) the bias disappears. We have also found that 
the Planck HFI experiment will be able to constraint the angle (3 with a precision of the order a{f3) ^ 0.1° at 95% 
and we have investigated the possibility to disentangle the systematic or cosmological origin of the rotation. However, 
the derived sensitivity on /? is one order of magnitude smaller than the precision achievable from direct, in-flight, 
measurements. It will be therefore mandatory either to calibrate the Planck polarimeters with better than 1° degree 
accuracy, either to include in future MCMC analysis /3 as an extra parameter. 

For the EPIC proposal, a value of /3 ~ 0.5° could heavily affect the cosmological parameters, if not properly included 
in the analysis, leading to a misleading claim of gravity waves detection at 3cr level. On the other hand the sensitivity 
on /3 will improve by one order of magnitude with ct(/3) ~ 0.01° at 95% c.l. in MCMC analysis. 
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